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1,    Introduction 

We  consider  in  this  paper  a  queue ing  process  that  was  encount- 
ered in  the  course  of  a  military  study.   Imagine  the  following 
situation;  A  missile  battery  is  defending  itself  against  a  train 
of  hostile  aircraft.   Each  attacking  plane  has  perfect  accuracy 
and  consequently  will  destroy  the  battery  if  it  gets  within  a 
certain  critical  range    .  V/e  thinly  of  the  file  of  aircraft  as  a 
queue  which  is  moving  toward  an  objective  0  and  of  the  battery  as  a 
server  who  attempts  to  perform  some  operation  on  each  element  in  the 
queue  before  that  element  reaches  0,  When  a  member  of  the  queue 
has  been  served  successfully,  the  server  retracts  instantly  to  the 
next  one  in  line.   If  an  element  in  the  queue  reaches  0  without 
having  previously  been  served,  the  server  is  permanently  and 
totally  disabled, 

A  similar  problem  was  formulated  by  B,  McMillan  and 
J,  Rlordan  [31 •   In  their  treatment  the  time  axis  is  regarded  as 
continuous,  and  the  service  time  is  assumed  to  be  an  exponentially 
distributed  random  variable.   Since,  in  the  aforementioned  example, 
each  la\inching  of  a  missile  renders  the  launching  facility  in- 
operable for  some  fixed  period  of  tirae,  independently  of  the 
success  or  failure  of  the  firing,  it  is  of  interest  to  study  the 
process  under  the  assumption  that  commencement  and  completion  of 
service  can  occur  only  at  definite,  regularly  spaced  time  points. 
Accordingly,  in  this  paper  time  is  treated  as  a  discrete  variable. 


71) 

This  idealization  is  made  to  simplify  the  model.   The  more 
realistic  case,  when  arrival  within  critical  range  results  in 
disability  of  the  battery  with  probability  a,  0  <  a,  <  1,  is 
treated  by  VJ,  Hirsch  and  H.  Kanish  in  "Some  Properties  of  a 
Transient  Queue"  (to  appear). 


*:i;^     i.'./l? 


2,   The  Ilathematical  Model 

Let  E- ,  Ep,  •••  denote,  respectively,  members  of  an  infinite 
queue,  i/here  the  subscript  Indicates  the  order  in  line,  Me   thinlc 
of  these  elements  as  being  located,  initially,  at  a  subset  of  the 
integer  lattice  points  on  the  positive  3t-axis,  the  precise  position 
of  each  element  being  stochastically  determined.   Specifically, 

let  4  be  a  random  variable  which  takes  the  value  k  with  probabil' ty 

co 

a,  ,  ^   a,  =  1,  and  suppose  that  at  time  t  =  0  (the  initial  instant) 

element  E^  is  at  the  point  £,    •   Tlie  distance  between  E.  and  E.,, 
1  ^      o  J      J+1 

is  a  random  variable  ^.  which,  for  all  j,  takes  the  value  k  viith 

J   .  .  .    .  ... 

probability  b,  ,  k  =  1,2^,,,  «   The  random  variables  I^i  ,  n=0,  1,,,A 
are  assumed  to  be  mutually  independent. 

The  queue  moves  with  unit  velocity  toward  the  origin  0,  its 
members  being  served  successively.   Service  on  E.  is  initiated 
instantaneously  upon  completion  of  service  on  E .  - ,   If  an  element 
arrives  at  0  unserviced,  the  server  is  disabled  (absorbed)  and  the 
process  terminates.   (No  interruption  occLirs  if  an  element  reaches 
0  and  simultaneously  service  on  it  is  completed).   The  tine  between 
commencement  and  completion  of  service  on  E.  is  a  randozn  variable 
X.   which,  for  all  j,  takes  the  value  k,  k  =  1,2,...,  l^^ith 
probability  C,  .   The  sequence  ^X   Ir  is  mutually  independent  and, 
moreover,  is  independent  of  -[^1  .   The  above  definitions  imply  that 
initiation  and  completion  of  service  can  occur  only  at  integer 
lattice  points. 
3»   Recurrence  Relation  for  the  Service  Distribution 

Let  V  denote  the  nuriiber  of  elements  served  before  termination 
of  the  ;jrocess.   Then  v  is  a  (possibly  improper)  random  variable 
and  we  set 


''■■y 
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and 


Pi^v  =  n|=  p^ 


P|v  =  nli^  =  kj  =  p{n,k) 


Clearly J 


CO 

(1)  Pn  =  m  P(n>l^)a 

^   k=l       ^ 


Let 

k 

Then,  since  v  =  0  if  and  only  if  E,  arrives  at  0  unserviced,  we 
have 

(2)  p(o,k)  =  qjj    . 

The   random  variable   v  takes  the  value  n  >  1   if  and  only   if 
service    on  E^    is   completed,   and  then  exactly  n  -  1  additional 
elements  are    serviced  prior   to  termination.      To  determine    the 
probability   of  tliis  event  vie  argue   conditionally,   assuming  that 
4     =  k  and  4-    =  m.      The   probability   of  completing  service    on  E,    at 
the   point    j,    0<   J<k-1,    is   simply. 


(4  =  k  -  j)  =  c,^..  . 


h 


If   this   occurs,    the    probability   that  exactly  n  -  1   additional 
elements  are    serviced  prior  to   the   absorption  of  the    server   is 
p(n  -  1,    j   +  ra).      Hence, 


k-1 

EI 
J=0 


)|^v  =  nU^  =  k,  c,^  =  mj  =  XZ  Cj^_jP{n  -  1,    j   +  m)      , 


and  therefore 


k-1     oo 
(3)        p(n,k)  =  !___  C      .YZ   P(n  -  1,  j  +  m)b 

j=0  ^  Jm=l 


The  rectirsion  relation  (3)  and  the  initial  condition  (2) 
Toniquely  determine  the  conditional  probabilities  p(n,k), 
n  =  0,1,,,,  •   The  absolute  probabilities  ^ pi  can  then  be 
computed  from  (!)• 

L|.,    Uniform  Spacing  and  Geometric  Service  Time 
Ij.,1  The  exact  solution.   Let  us  now  assume  that  the  time  spent 
by  the  server  on  a  given  queue  element  is  determined  by  a  sequence 
of  independent  attempts  at  service,  each  taking  unit  time  and 
having  probability  p  (o<  p  <  1)  of  success «   The  service  time  V. 
is  then  geometrically  distributed,  lee.. 


c^  =  P.'^-^ 


where  q  =  1  -  p.   Suppose  further  that  the  distance  4 •  between 

J 

consecutive  elements  in  the  queue  has  the  degenerate  distribution 


:.'ri 


p{53  =  d]  =  l 


where   d   is   a  positive    integer. 

The   initial  condition   (2)  becomes 


ik)  P(0,k)   =  q^  =  q^   , 


and  the   reciorslon  relation   (3)   takes  the    form 

(5)  p(n,k)   =  Si:Pq^''"'"^P(n-l,    j+d)      . 

0=0 

The  v_nique  solution  of  {$)  xander  the  constraint  (Ij.)  is  given  in 
the   following  theorem: 

Theorem  1«      If   the   elements   in   the  queue    are    spaced  at  equal 

k-1 
intervals   d,   and  the    service   time   has   the   distribution  C,    =  pq        , 

then 

(f.)  r^fv.   1^^   _     k(k+nd-l)l   •       n„k+n(d-l) 

(6)  P^^'^^)   =nilk+n(d-l)jj    P  ^  V        , 

Proof,  The  proof  is  inductive.  It  is  evident  from  {I|.)  and  (5) 
that  (6)  holds  for  n  =  1<>  Suppose  it  is  true  for  n  =  r»  Then, 
defining 

we  have 


i   '-%  -'  <'H 


r      ,,^.-1.'"' 


1.  '  ■  .  fi  /  . 


r+l  k+(r+l)d-l)  k-1  ,^. 

p(r+l,k)  =£ ^, n:(J+<a)C0+(r+l)d-l]^^-^^ 

rMa^k+(r+l)(d-l)  k^ 


Using  the  identity 


.-^^ —  r  +  1 

a=j:n 


we  find  that 


-  d[k+(r+l)d-ii^^V 

r+1  k+(r+l)(d-l)  ,    \ 

==^ Vin k[k+(r-:-l)d-l]^^)   , 


which  completes  the  induction. 

I|.«2  Asymptotic  behavior.   The  limiting  behavior  of  p(n,k)  as 
n  — >  oo  is  obtained  easily  by  using  the  Sterling  approxiraation 
and  the  relation 

(n+m)!  ^-^  nin^       • 

For  d  =  1, 


/      ,  X       ^  ■     „k-l  n 


while   for  d  >  1 


cp^ 


p(n,k)  — '--!?72      * 


where 


C   =  k  '  ^ 


■■'^"^J  v^dTdTry  ' 


and 


6  =  i;a. ^L^    . 

^   (d.l)'^-l  • 

The  function  p  =  P(p),  which  takes  values  in  the  interval  0  <  P  <  1 
when  0  <  p  <  1,  is  monotone  increasing  on  0  <  P  <  l/d  and  monotone 
decreasing  on  l/d  <  p  <  Ic   Clearly, 


P(|)  =  1  . 


Consequently,    if  p  <  p"   <  l/d,   vie   have    (in  an  obvious  notation) 

p(n,k)    =    o(p''(n,k)) 
as  n  — >oo.      If  -T  <  p  <  p   , 


p"(n,k)  =  o(p(n,k) )      . 


8 

On  intuitive   grounds  vje  might  have   anticipated  that  for   sxifficiently 
large   values   of  n, 

p(n,k)   <  p'(n,k)      • 

Paradoxically,  however,  the  opposite  inequality  holds  in  the  ran^e 
-T  <  p  <  1,   The  apparent  paradox  is  resolved  in  the  next  sectio.ia 
where  it  is  shown  that  there  is  a  positive  probability  of  servicing 
infinitely  many  eleinonts  if  and  only  if  p  >  l/d,  so  that  for  values 
of  p  in  this  range  some  of  the  probability  mass  is  shifted  out  to 
infinity, 

I|.,3  The  probability  of  vtltimate  absorption.   In  this  section  we 
study  the  behavior  of  the  quantity 


-22- 
f  =  1 p(n,l:)   , 

n=0 


which  can  be  interpreted  as  the  probability  that  the  server  is 
disabled  ultimately,  l.eo,  after  serving  some  finite  number  of 
elements.   (The  difference  1-f  is  the  probability  that  all  members 
of  the  queue  are  served,  or  equivalently,  that  the  server's  life- 
span is  infinite.)  For  this  purpose  it  is  convenient  to  introduce 
the  generating  function 

oo        . 

F(:cjk)  =  5   p(n,k)x^  , 

n=0 

restricted  to  tlie   domain  0  <  x  <  1.      The   properties   of   P  of   inter- 
est  to  us    can  be    obtained  expeditiously  by  using  the   following 
characterization: 


'.    i        *•>. 


I'HrlJ 


Theorem  2.   If  p  ^  ^,  to  each  real  nuiaber  x  in  the  closed  interval 
0  <  X  <  1  there  corresponds  exactly  one  value  z  =  z(x)  in  the  disc 


N-1 1  <  dti 


1  (2) 


such  that 


(7)  z  =  1  +  xpq^'^z^  . 


The  functions  P(x,k)  and  z(x)  are  connected  by  the  relation 


(8)  P{x,k)  =  [qz(x)]^^  • 


If  p  =  ~>,   the    same   result  holds   in  the   half-open  interval 
0  <  X     <  1. 


Proofs      If  d  =  1,    equation   (7)   has    the   unique    solution 


z(x)  =   ^ 


1-xp  * 
which  by  bhe  binomial  expansion,  leads  to 

n=l      '^* 

Comparing  this  equation  with  (6),  we  see  that  for  d  =  1 

__ ^ 

If  d  =  1,    the   disc  may  be   talcen  to  be   the   infinite    z-planeo 


x"      • 


10 


P(x,k)    =    [qz(x)]^  =  U^    J 


Suppose   now  that   d  >  1.      Let   G  denote   the    open  disc 


jz.ll   <^l-«     , 


d-; 


and  let  f  be  its  boundary.   The  functions  z  and  —r-  are  anaj.-j^tlj 

—  z 

in  Gf   the  closure  of  G,  and  it  is  easy  to  verify  that 


inf. 

ze  r 


_d" 


d-l 


Applying  the    theorem  of  Burmann- Lagrange  we   conclude   that  to  each 
point  w   in  the   disc 


lwl<i^ 


d-l 


there  corresponds  a  unique  point 


in  G  such  that 


z  =  t})(w) 


z  =  1  +  v;z  ■  , 


and  the    function 


Z^    =    4'^'(V7) 


11 


has  the   series  representation 

1,                        oo_  k(k+nd-l)^^-^j.n     ' 
=  4>   (w)   =  1   +  ^^ --j ^         w 


-k  _   iki         ^  


It  folloivs  from  the   inequality  of  the  arithmetic   and 
geometric  means   that   if  either  pj^-rorx^l, 

d^ 

Hence  we  may  replace  w  by  xpq  ~      in  the  Burmann-Lagrange 
exnansion.      Setting 

(Kxpq"^-^)    =  z(x) 


and  multiplying  both  sides  by  q  we  find  that 

,     1    00   n„k+n(d-l),  z,^  ^  T\(n-1) 

(9)   [qz(x)]^  =  ql^  +JZ^-2 ll(|pd^lj ^n  ^ 

n=l 


where   z(x)    is   the   unique   root   in  G  of  the   equation 

z   =  1   +  xpq^"-'-z^     . 

We   observe   that   the   right-hand  side   of   (9)    is    P(x,k),  which 
completes   the    proof, 

Stippose  now  that   p  =  ■?•      Then  the   function  z(x)   appearing 
in  the   preceding  theorem  is  defined   only  on  the    half -open 
interval   0  <  x  <  1,     However,    its   definition  can  be  extended  to 


r 


I- 


x)-|> 


'•-   ri-tocf  ■Qcilv 


i^i 


,      "s        pqx  +   I    "-   s 


;    ^O    3: J. I  a 


3il&  no  rXno  b&al-'ieb  zt  rv^'ioe'r. 
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the  closed  Interval  so  that  relations  (7)  and  (8)  are  preserved  by 
setting 


(10)  z(l)  =  lira  _  z(x)   , 

z — >l' 


The  existence  of  this  limit  is  guaranteed  by  the  continuity  and 
monotonlcity  of  z«  From  Abel's  theorem  we  conclude  that 


P(l,k)  =  [qz(l)]^  , 


Moreover,  letting  x  — >  1"  in  (7)  it  is  clear  that  z  =  z(l) 
satisfies  the  equation 

Mt  \           ,  -  1  +  v.n^-l.d  -  l+(d-l)'^"-'-  ,d 
(11)  2  =  1  +  pq   z   =  ■  ^ Z    • 

d^ 

We  cannot  assert,  however,  that  the  value  z(l)  defined  by  (10)  lies 
in  the  disc 


z-ll  <^   . 


Indeed  it  is  easy  to  see  that  (11)  has  the  unique  real  root 


^  =^  -^dTi   » 


which  lies   on  the   bo\indary   of  the   disc,    since   the   function 


Z 
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increases  rtionotonically  in  the  interval  0  <  z  <  1  +  XTt"*  attains 

/  ,  ,  \d-l  ■> 

the  value   "^  "'•         when  z  =  1  +  -r-yf   ^^'^'^-   thereafter  decreases 

monotonically* 

We  nov;  use  the  representation  of  P  given  above  to  derive  a 

necessary  and  sufficient  condition  that  f  =  1: 


Theorem  3.«     The  probability  f  of  ultimate  termination  of  service 

has  the  value  1  if  and  only  if  d  <  — ,  i.e.,  the  distance  between  ad- 

jacent  elements  in  the  queue  does  not  exceed  the  expected  service 

time* 

Proof.   If  1  <  d  <  — ,  the  point  —  lies  in  the  disc  |z  -  ll  <  XTT* 

However,  z  =  —  is  a  root  of  (7)  when  x  =  1.   Since  there  is  but 

one  I'oot  in  the  disc,  we  have 


z(l)  =  l/q 


and  therefore 


f  =  P(l,k)  =  Cqz(l)]^  =  1 


If  d  =  — ,  equation  (11)  has  the  unique  real  root  z(l)  =  —  and, 
once  again, 


(3) 

'-''After  this  theorem  was  obtained  it  was  shown  by  W,  Hirsch  and 
H.  Hanish  in  "Some  Properties  of  a  Transient  ^.ueue"  (to  appear) 
that  it  is  a  special  case  of  a  much  more  general  result;   In 
a  wide  class  of  queueing  processes  with  an  absorbing  state 
ultimate  termination  of  service  is  almost  certain  if  and  only 
if  the  expected  distance  between  adjacent  elements  in  the 
queue  is  not  larger  than  the  expected  service  time. 


•!.■  T 


-■■i 


'-£?-^ 


«t  V  i-X^-' 


o  '-%  »*  >.s  -■ 


-jCiV 


er--M 


■'■l^a 


^''{Ql 


Ik 


f  =  1 


To  complete  the  proof  we  note  that  if  d  >  — ,  then  q  <  -^*  and 


f  =  [qz(l)]^  <  [-^(1  +  -^)]k  =  1. 


When  f  <  1,  it  is  possible  to  estimate  its  magnitude.  We 
do  this  by  approximating  z(l)^ 


Ay 


(d-1 

,d 


\d--l 


pq 


d 
.d-1 


(z^',y(z^)) 


Since 


z(l)  = 


P(l,k)  1^  _ 


L  q 


^z 


^  \r   ,22-. 

1+q   EZ  P(n,k) 

n=l 


1 

IE 


>  1  , 


a  first  crude  estimate  Is  given  by  the  interval 


1  <  z(l)  <  1  +  -^   . 
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_  z~l 


Now  y(z)   =  — y-  is   concave    In  this   interval,   and  consequently  the 

z*^  11 

chord  C^(z),    connecting   (1,0)    to    (1   +  "T^*  yd   ■*''drr^^»   lies 

completely  below  y(z)    in  the    interval.      By   solving  the  equation 


Cn  (  z  )    =  pq 


d-1 


for  z  we  obtain  an  improved  upper  boiind  for  z(l)«   Let  us  call 
this  momber  z,  ,   If  we  now  use  the  chord  Cp(z),  connecting  (1,0) 
to  (z-,y(z, )),  and  repeat  the  described  process,  we  obtain  a 
further  refineraent  "zp.   This  iteration  scheme  yields  a  decreasing 
sequenced  "z  V  converging  to  z(l).   The  values  of  s  "z  (  can  be 
obtained  recursively  from  the  relation 


"^  =  1  +  pq^"^(z^_3_)^,   n  =  1,2,..., 


with 


Z^  =  1  +  -~r         . 

o       d-1 


pq 


d-1 


>^  od"  nniatevaoo 


,-....%.t  ^  a       ,^(        "■•"■^pq  +  X  ==  jj" 


(':)V 


<  - 


f 


/ 


/ 


/ 


r^ 


t  ! 


\ 


\ 
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We  improve  the  lower  bound  for  z(l)  by  tangential  approxim- 
ation to  the  curve  y(z),  A  first  improvement  is  obtained  by 
constructing  T, ( a ) ,  the  tangent  line  to  y(z)  through  the  point  (1,0), 
and  solving  the  equation 


T^(z)  =  pq^-^ 

for  z*  Let  z.  denote  the  solution.  For  further  refinement  we 
construct  Tp(z),  the  tangent  line  through  iz^  tjiz^) ) ,  and  find 
Zp   from  the  equation 


T2(z)  =  pq^-^   . 


This  iteration  process  yields  an  increasing  sequence  f ^  ^  converging 
to  z(l)»   It  is  easy  to  verify  that 

.  _  P^    "^n^l^-^-n-l-^^^^^^^  -^^-n-l     „  -  .  o 

^  z-%  -  dz-%-l(z   ,-1) 

— n-1    —n-l  — n-1 


with 


^0  =  1  • 


ij-.I).  Moments  of  the  service  distribution.   By  successive  differ- 
entiation of  the  equations 


and 


P(X,1C)  =  [qa(K>]^ 


z(x)  =1  +  xpq'^'^z^Cx) 
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the  moments   of   v,   the  nujnber   of  elements   served  before   termination, 
can  be   calciilated.      The   first  two  of   these   are   given  in  the 
following  theorem: 


Theorem  k"     The  expected  number  of  elements   serviced,   E[v],    is 
finite    if  and  only   if  d  <  --.      In  this   case 

and 

(1-dp)-^ 
Proof.  Using  the  algebraic  equation  satisfied  by  z(x)  we  find  that 

Zi\.X; -rr^s        -r^  • 

l-xpq"^  -"dz"^  -^(x) 


Since 


Pt(x,k)  =  kq^S^^"^(x)z«(x)   , 


we  have 


(12)  P.(x,k)  =  ^ ,  ."   ,  ,^-^^ 

l-xpq'^'-^dz^-  -^(x) 


It  was  shown  in  the  proof  of  Theorem  3  that  z(l)  =  l/q  if  d  <  —  . 
Setting  X  =  1  in  (12)  it  follows  that 


E[v]  =  p'(l,k)  =^     . 


18 

If  d  >  — ,  Theorem  3  shows  tliat  Pjv=  oo?  >0,  and  hence  v  has 
infinite  expectation.   Thus,  the  expected  value  is  finite  if  and 
only  if  d  <  ^. 

The  variance  is  calculated  by  a  second  differentiation.  For 

P 


d  <  —  we  find  that 


P"(l,k)  =-i^(k  4-d  ^^) 
(1-dp)^  ^  "^P 


Hence, 


<tHv]   =  p"(l,k)  +  p'(l,k)  -  [p»(l,k)]2  =      ^^3  . 

(1-dp)-' 


5 ,   C-eonetric  Soaping  and  Geor.etric.  Service  time. 

5ol  T'-ae  Generatin"  Function 

In  this  section,  as  in  the  previous  one,  we  assume  the 
service  tirae  to  be  geometrically  distributed  with  parameter 

p(0<p<l),  so  that 


0^  =   pq^"^  ,   k  =  1,2,..,   ♦ 


The  distance  betvjeen  adjacent  elements,  however,  is  no  longer  taken 
to  be  uniform.  V.'e  assurae,  rather,  that  it  is  determined  by  a  .  . 
succession  of  Bernoulli  trials,  each  having  probability  "p  (0<  p"  <  1) 
of  success.   The  distance  takes  the  value  k  if  the  first  success 
occurs  at  the  kth  trial.  Hence, 


b^  =  pq  ■   ,   k  =  1,2,,.,   , 


'■^ 
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Although  in  principle  the  probabilities  '?p(n,k)j'  can  be 
calculated  recursively  from  (2)  and  (3)>  it  is  difficult  to 
obtain  a  general  expression  in  this  xijay.   On  the  other  hand,  the 
explicit  form  of  the  generating  fuaiction  P(x,k)  can  be  obtained 
with  comparative  ease.   Using  (2)  and  (3)  we  find  that  P(x,lc) 
satisfies  the  equation: 

,       k-1  ,  .  ,  oo  - 

(13)       P(x,k)  =  q^  +  ppx  tH   q^"^   IZ  P(x, j+ra)q^--^  . 

j=0        m=l 

Setting  X  =  0,  ^^re  have 


k 


(14)  P(0,k)  =  q 


Moreover,  successive  differentiation  determines  the  values 


d"P 

dx" 


at  X  =  0,   Hence  (13)  has  only  one  solution  which  is  analytic  in 
the  \anit  disc. 

If  we  substitute  a  trial  solution  of  the  form 


P(x,k)  =  U^(x)  , 


we   are   led  to  a  g_uadratic   equation  for  U; 


qU^   +    (ppx   -l-^q)U+q=0      . 


The  root  satisfying  the  condition  U(0)  =  q  is 


I  ,i 
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(15)  U    (X)    =   1    +   qn    -    PPX    ~    fd-qo)^  ■    2pB(l+qq)x    +    (pp)^X^3^/'^    . 

Since   this   function  is  analytic   in  the   unit   disc   and   satisfies    (13 )# 
it   follov;s   that 


(16)  P(x,k)   =  ul^(x)      . 


Putting  X  =  1  in  (l5)  we  see  that 

U3_(i)  =^  ^^,:  1^  -^1    , 

which,  together  with  (16),  establishes  the  following  theorem: 


Theorem  5»^   The  probability  f  =  P(l,k)  of  ultimate  termination 
of  service  has  the  value  1  if  and  only  if'q<q,  iee«,  the  expected 
distance  between  adjacent  elements  does  not  exceed  the  expected 
service  time.   If  "q  >  q^ 


^-qr'  . 


The  moments  of  V  are  readily  calculated  from  (l5)  and  (l6). 
We  give  the  first  two  of  these  below: 

Theorem  6.   The  expected  number  of  elements  served  is  finite  if 
and  only  if  "p  >  p»   In  this  case 

^  This  is  another  instance  of  the  general  situation  described  in 
the  footnote  on  p,  13 • 


l£Ji.iJ 


P  *s  P  'il 
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E[v]    =   P«(l,k)    =  ^     , 

P-P 


and 


cr^Cv]    =   P"(l,k)    +   P'(l,k)    -    [P'(l,k)]^   =  — ^^(p^q    +   p^q) 

(p-p) 

5«Z  Asymptotic  behavior >  The  behavior  of  the  probabilities 
fp(n,k)j  for  large  values  of  n  Is  described  In  the  follox^/lng 
theorem: 


Theorem  7,     As  n  — >co, 

p(n,k)^^2^^^      , 


where  C  Is  a  constant  (depending  on  k),  and 


J.  =  Ppd-'-^iq-S  y/qq) 
(1-qq)^ 


Proof.   In  the  relation 

P(5  ^O  -A  +  qq  -  PPS  -  [(l-aq)^^2pp(l4qq)s  +  (py)^S^]^^n^ 
L  ^  ^  J 

make  the  substitutions 

i-q<t 
(17)  x=:^s   , 

l-qq 
y(x)  =  (2qr)^(l-^q)-^^P(o,k)   . 


m 


It  is  Interesting  to  note  the  similarity  between  this  result  ^nd 
the  corresponding  one  in  the  case  of  uniform  spacing. 


Then 


•(x)    =  /a   -  X   -    VTI2H+?  j^      , 


Taking  the   logarithmic   derivative. 


2Z 


(18) 


I    = 


Ja. 


Vi-2i\x+x^ 


and 


y"   =  k 


Jll 


2T372 


N/TIil^?        d-SAxfx'^) 


Hence  y  satisfies  the  differential  equation 


(19) 


(1  -  2Ax  +  x^)y"  =  k^y  -  (x  -  A)y«   , 


If  we  express  y  in  the  form 


(20) 


,   GO   A^ 

y=(A.l)k^^.n   ^ 


by  equating  coefficients  in  (19)  v;e  find  that 


(21) 


^n+2  =  ^^2n+l)A        +   (k^-n^)A^,        n  =   0,1,...      . 


n' 


Since 


y(0)    =    [2q(l-qq)]-kx-(0,k)    =    [  2q  (l-qq  )  ]  "^qk   ^  f_M)k   ^    (^.^^k 

\i-qq/ 
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it   follows  from   (20)    that 


(22)  ^o  "  ^      • 


Similarly,  ue  see  from  (18)  that 


y'(0)  =  1.-5^(0)  =  k(A-l)^   , 


and  hence  from  (20) 


(23)  A^  =  k   . 


We  now  put 

C^(n+k-l.)J 
(2^)  A^  =  ^-.^^^     . 

Then  (21)  becomes,  after  simplification, 

(25)       (n+k+l)C^^2  =^(2n+l)C^^^  +  (k-n)C^   , 

and  the  initial  conditions  (22)  and  (23)  transform,  respectively, 
into 


and 


^0=1 


C^  =  1   . 


To  solve  (25)  we  use  the  familiar  method  of  Laplace.   Set 


2h 


b 


'n+1 


C   t''"^[tix(t)]dt      , 


^n+2  =\  t'"^[t^ix(t)]clt   , 


nCn  =  t   \i{t) 


b       b 


-  C  t''"^[t[i'(t)]dt      , 


a       a 


and 


(n+l)C^^3_  =  t^'*"^tx(t) 


(n+2)C^^2  =  t^'"^l^(t) 


b       b 


-V  t'^'^Lt^iJ.'Ctjldt     , 


a       a 


b    .    b 


-  [  t^'^[t^ix«(t)]dt      , 


a       a 


where   a,   b,   and   [i{t)   are   to  be   determined,     VJith  these   substitutions 

(25)  takes   the   form 

(26)  rt^~-^Yt^-2At^+t)ia'(t)    -    [(k-l)t^+At-k]tL(t)ldt 


-  t^(t^-2At+l)ii(t) 


=   0      . 


To  satisfy  (26)  it  suffices  to  determine  p-(t)  so  that  the  coefficient 
of  t^~  is  0  and  then  to  choose  a  and  b  so  that  the  integrated  term 
is  0,  Accordingly  ^  satisfies  the  differential  equation 


4---.    ..X.-a. 


l(u  )'>Lj  -J  ]' 


iij 
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(27) 


ILL  =   (k-l)t^   +  At-k  ^ 

^       t^-2/a^+t 


^'  +   n.   In    •'2A+2t 


Solving  (27)   xre    obtain 


ix(t)    = 


?  V  1 


The   equation  for  the    integrated  terra. 


(28) 


t'^(t^-2At+l)^'^2 


=    0 


is  satisfied  if  a  and  b  are  roots  of  the  quadratic  equation 

t^  -  2At  +  1  =  0   , 

If  n  >  k,   we   can  also  choose   a   =  0,      Thus,   for  n  >  k,    {28}   has 
distinct   independent   solutions 


a 


/  'p 

=   0,  b   =  A    -    V  A   -1 


and 


a  =  A  -    \/i?-l    ,   b  =  A   +  Va^-1      , 


Hence,    if  n  >  k. 


(29) 


^n  =  h 


^  t^-k-l(3__2At+t2)'  ^dt 


k-4 


11, 


^  t^-k-l[.(l_2At+t2)]    "2^t      , 
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where 

r^  =  A  -    l/?-l        , 

r^  =  A   +  i/a^-1        , 

and  K, ,  Kp  are  constants  independent  of  n,  (The  negative  si^n  is 
used  in  the  second  integral  to  make  the  terra  in  brackets  positive 
in  the   interval   of  integration. )      Noting  that 

f  2  k  i  ^1  k  J= 

\    t^-k-lfi.2At+t^|      ^dt  =(    t""^'-^(l-2At+t2)      ht 


0 


(    t--^-l(l-^ 


^2  _   ,.  ,  o     k4 


+ 
r. 


C    t^-l^-^[.(i.2At+t^)]    "^dt    , 


we   re-write    (29)    in  the   form 


(30)  C^  =  l^{    t^-^"^|l-2At+t2l  '   ^dt   +  t^  [    t"-^-^|l-2At+t2|^'"^dt, 

0  'o 

where  £-,  and  Oo  are  constants  independent  of  n. 

The  asymptotic  behavior  of  these  integrals  as  n  — >oo  is  re- 
vealed by  making  the  transformation 


t  =  r^e"""  ,   i  =  1,2  ♦ 


Then  (30)  becomes 


.'  ■'  f  1 


•,t--av 
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1    00 


n 


J-      OO  V     J;  l'_i 


0 


1     OO 

i 

0 


k-1 


VrX«''""""li-^'"l      1^1-2^'"! 


-PpS      I        au» 


Let 


k4 


a 


F^(u)    = 


U-e-^l      2|r     r^e-^l      ^ 


and 


1     GO 


f^(n)   =  vl  2  Je-^^-^)V(u)du        , 


0 


where    1  =  1,2,   and   j   =   j(l)   = 


_  A 


2   if   i  =  1  . 


•      Since 


k-4 


1   If   1  =  2 


k-i 

-U|i      2 


F^Cuj-^u        Ir.-f^e'    f 


as  u  — >  0,    it  follows   from  the   fojrailiar  Abelian  theorem  for 
Laplace   transforms  that   as  n  — >oo 

.    1 


k-i 


1 


k-4 


fi(n) 


r.   ^[r.-r.|      ^fCk^)  r.   ^Ir.-r.T  ^Hk^) 


(n-k) 


kH 


n 


H 


Moreover,    since  r,    <  Vp, 


f3_(n)   =    0(^2^^)) 


W) 


Cf   [?],    P  181,   Theorem  1,   Corollary  la. 


oo 


^~i.l 


.:>    ^-l 


> 


'    -t  cJ'i 


:s. 


^S*^ 


•      K{n)r^'i)o  '   («)x'i 
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Hence, 

(31)         C^^-^^ ^  ^  ., ^— •   . 

n 

Transforming  back  to  the  original  variables  (see  (17),  (20),  (21|), 
and  ( 31 )}.,  we  have 

/   1  \    /^-5l5'Nk/A    .,  sk/  pp  ^n   (n+k-l)l^  Cr^ 

p(n,k)  =  {-^)    (A  -  1)  (3^)   ^^(^,i)j-C^-  -372   , 


where 


=  !i!?  =  £iLLt2Sz2iMi  ^"^^ 


1-qq       (l-qq) 


This  completes  the  proof, 

TtT 


It  is  easy  to  verify  that  r  <  1.   Thus, 


r  pp   r2PP      ]^-  ^  q+q,  ,^,  ^-i 

1-qq   1-qq    '^  1-qq    1-qq     '^  ^ 


<  r2[A    -    \/qq(A+l)]    =   r2[A    -    jf3(A+l)] 


=  r2(A   -VA^-l)   =  r^r2  =  1, 


'.•- 
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6,   Remarks  on  the  Continuous  Case 

It  is  interesting  to  note  that  results  for  the  case  when  time 
is  treated  as  a  continuous  variable  can  be  obtained  by  an 
appropriate  passage  to  the  llrait  in  the  corresponding  discrete 
forrnuli.   To  this  end  make  a  change  of  scale  so  that  the  discreto 
process  takes  place  on  the  lattice  points  \  mh,  m  =  0,1,...  i   • 
Suppose  the  first  element  in  the  queue  is  at  a  distance  kh  from 
the  absorbing  barrier  0,  and  the  distance  betX'Jeen  consecutive 
elements  is  dh.   Then,  as  the  analysis  in  §l|.l  shows. 


p(„,k)  =  Mkina^lji"li!  pn^k-^ntd-l)   , 


Now  let  a,  T,  and  e  be  positive  constants,  and,  as  in  the  passage 
from  Bernoulli  trials  to  the  Poisson  distribution,  let  h  — >   0, 
p  — >  0,   k  — >oo,  and  d  — >co  so  that  kh -v_^T,  dh  ^-^e  and  p^--^ah. 
Then 


pU,kJ  '"^ — i 


which  agrees  with  the  result  in  [3]j  where  the  service  time  is 
exponentially  distributed,  the  initial  element  is  at  a  distance  T 
from  0,  and  the  spacing  is  uniform  at  intervals  e. 

If,  in  the  discrete  model,  the  spacing  is  geometric,  let 
'p',w  ph,  where  p  >  0,   Then,  for  the  generating  function  P(x,k)  we 
obtain,  in  agreement  x^jith  [3],  the  asymptotic  formula 


'■>    r  — . 
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^  T  . 


Pormuli  for  the  moments   of  the   niAmber  of  elements   served 
before    termination  can  be    obtained   in  the   same  way. 
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